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Abstract 

The aim of the article is to discuss the S'-matrix interpretation of perturbation 
theory for the Wigner functions generating functional at a finite temperature. For 
sake of definiteness, fruitful from pedagogical point of view, the concrete problem 
from particle physics of high-temperature initial states dissipation into cold one is 
considered from experimental and theoretical points of view. The temperature is 
introduced in the theory by typical for the micro canonical description way. The 
perturbation theory contains two- temperature (of initial and final states) Green 
functions. Two possible boundary conditions are considered. One of them is usual 
in a field theory vacuum boundary condition. Corresponding generating functional 
of Wigner functions can be used in the particle physics. Another type of the bound- 
ary condition assumes that the system under consideration is in environment of the 
black-body radiation. This leads to the usual in statistics Kubo-Martin-Schwinger 
boundary condition at the equilibrium (one-temperature) limit. The comparison 
of the S'-matrix approach with Schwinger-Keldysh real-time finite-temperature field 
theory and with nonstationary statistical operator approach of Zubarev are consid- 
ered. The range of applicability of the finite-temperature description of dissipation 
processes is shown. 



1 Introduction 



At the very beginning of this century couple P. and T.Ehrenfest had offered 
a model to visualize Boltzman's interpretation of irreversibility phenomena 
in statistics. The model is extremely simple and fruitful |Ij]. It considers 
two boxes with 2N numerated balls. Choosing number / = 1,2, ...,2iV 
randomly one must take the ball with label I from one box and put it to 
another one. Starting from the highly 'nonequilibrium' state with all balls 
in one box it is seen tendency to equalization of balls number in the boxes. 
So, there is seen irreversible^ flow toward preferable (equilibrium) state. 
One can hope QXQ that this model reflects a physical reality of nonequilib- 
rium processes with initial state very far from equilibrium. A theory of such 
processes with (nonequilibrium) flow toward a state with maximal entropy 
should be sufficiently simple to give definite theoretical predictions. 

In order to do the consideration less formal we will be connected with 
concrete physical problem. For instance, the particles creation processes 
are good laboratory for investigation of relativistic nonequilibrium pro- 
cesses general properties. Indeed, considering the multiplicity n as the 
characteristics of final state entropy we can choose the asymptotically large 
n » n(s), where mean multiplicity n(s) naturally defines the scale of n. 
Then one can expect, noting above mentioned general property of the 
nonequilibrium flow, that the theory of processes with practically total 
dissipation of initial-state kinetic energy into particles masses should be 
extremely simple. By this reason it is natural to start consideration from 
region n » n(s). We would construct the theory permanently taking into 
account just this condition. 

The theory of dissipative processes have general significance from ther- 
modynamical point of view and we would concentrate attention on this 
important problem. There is also practical side of considered problem. At 
n >> n(s) the cross sections a n (s) falls down rapidly and are too small 
(< nb). Noting also a problem of triggering such rear final state the ex- 
perimenters must have enough arguments to examine them. The main 
arguments are as follows: at n >> n(s) we have unique chance (i) to ex- 
amine 

11 What never? No never! What never? Well, hardly ever.' || 
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- pure (practically without admixture of hadrons), 

- cold (it is a best condition for investigation of collective phenomena in a 
system) , 

- dense (in this case the QCD interaction constant a s is small) 

quarks plasma (CQGP) and (ii) realize experimentally decay of very hot 
(at high energies) initial state in the 'inflational' regime, with 'freezed' 
nonperturbative degrees of freedom of hadrons system. 

It is known from hadrons high-energy experiments that the cross sec- 
tions o n have a maximum at n ~ n{s), where 1 << n(s) << n max and 
n-max = y/s/mh is the maximally available multiplicity at given energy 
■s/s (rrih is the hadrons characteristic mass). This testify to the state- 
ment that in hadron processes the nonequilibrium flow is not equal to zero 
(n(s) >> 1), but the mostly probable processes did not lead to the state 
with maximal entropy (n(s) << n max ). (The early models was based on 
the assumption that the final state of inelastic hadron processes has max- 
imal entropy n(s) ~ n max g.) 

The preferable at n ~ ft processes are indebted for excitation of hadrons 
nonperturbative degrees of freedom described by creation of hadrons con- 
stituents from vacuum: the kinetic motion of partons leads to increasing, 
because of confinement phenomenon, polarization of vacuum and to its 
instability concerning quarks creation ||]. In other words, there is a long- 
range correlation among hadrons constituents. Under this special correla- 
tions the conservation laws constraints was implied. They are important 
in dynamics since each conservation law decrease number of the dynamical 
degrees of freedom at least on one unite, i.e. it has nonperturbative effect 
(this must explain why n(s) << n max ). Moreover, in so-called integrable 
systems each independent integral of motion (in involution) reduces num- 
ber of degrees of freedom on two units. In result there is not stochastization 
in such systems H, i.e. the nonequilibrium flow is equal to zero. But it 
will be argued that at the very high multiplicities this effect is negligible. 
So, if 

n(s) «n< n max 

we will see that the particles creation processes are close to Markovian in 
accordance with Boltzman's idea. The reason of this phenomena is the 
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more fast falling down of soft chanel of hadrons creation comparing with 
hard channels in asymptotics over n. 

Rejecting nonperturbative effects creation of the high-multiplicity final 
state can be described using standard methods of QCD. We will show 
dominance of processes with minimal number of QCD jets in the high- 
multiplicities region. This means that the high-multiplicity processes are 
stationary Markovian^. This result is in agreement with Boltzman's general 
idea concerning nonequilibrium flows. 

So, the high-multiplicity processes are 'unshadowed' by nonperturba- 
tive and complicated perturbation effects. This will allow to investigate 
not only new state of the pure colored plasma but also the structure of 
fundamental Lagrangian. This conclusions are not evident and we start 
consideration from brief review of arguments. 

It must be noted that the experimental investigation of high-multiplicity 
processes in deep asymptotics over n seems unreal. But considering mod- 
erate n >> n we can not be sure that the final state is equilibrium. In- 
vestigation of fractal dimensions in the multiparticle hadron processes at 
high energies shows presence of considerable fluctuations Q5j . This leads to 
necessity to have the theory of dissipation processes with nonequilibrium 
final state. 

There is also another side of the problem. Today understanding of 
hadron processes is far from ability to give any quantitative prediction. 
The above announced prediction concerning absence of nonperturbative 
contributions into hadron processes 'work' in the deep asymptotics over 
n only. So, at moderate n >> n we can not be sure that they do not 
have important influence. That is why we will concentrate our attention 
in this paper on the searching of economic (thermodynamical) description 
of the dissipative processes, trying to find the connections of our iS-matrix 
approach with other ones. It is important to note that the offered for- 
malism allows to separate the dynamical side of question from the pure 
descriptional one (see also concluding Section). 

This central problem of formalism can be solved noting that our dissipa- 
tive problem contains element of dynamics since it crucially depends from 

2 The vertexes renormalization takes into account the time-reversed fluctuations in the nonequilibrium 
flow. 
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boundary condition. Therefore, we adopt ^-matrix formalism which is nat- 
ural for dissipative systems time evolution description. For this purpose 
the amplitudes 

< (P)m\(q)n >= On,m(jPl:P2: —,Pm', Ql, Q2, —,q n ) 

of the m- into n-particles transition will be introduced. (The in- and 
out-states must be composed from mass-shell particles 0.) Moreover, 
to incorporate the boundary condition n >> m we should calculate the 
probability integrating over particles momenta: 

r(P;n.m) ~ J |a n ,m| 2 = J < (p)m\(q)n >< (q)n\(p)m > 

since the amplitude a nm (pi,p2, ■■■,Pm'i <7i> q2i ■■■■> qn) is the function of too 
many variables. This standard method of particles physics practically 
solves our problem. 

Nevertheless it is desirable to use the thermodynamical language as the 
mostly economic one, i.e. the formalism which uses minimal number of 
parameters (temperature, chemical potential, etc.) for description of the 
system. 

The field-theoretical description of statistical systems at a finite tem- 
perature is usually based on the formal analogy between imaginary time 
and inverse temperature (3 ({3 = j8j]. This approach is fruitful [|9j for 

description of the static properties of a system, but it demands a compli- 
cated mathematical apparatus for the analytic continuation to the real time 
TUf , if we want to clear up dynamical aspects. The first important quan- 



titative attempt to build the real-time finite-temperature field theory [JTT| 
discover a problem of the pinch-singularities. Further investigation of the 
theory have allowed to demonstrate the cancellation mechanism of these 
unphysical singularities \\F2§ . This attained by doubling of the degrees of 



freedom |T3|, [TJj] : the Green functions of the theory represent 2x2 matrix. 
It surely makes the theory more complicated, but the operator formalism 
of the thermo-field dynamics |T5j shows the unavoidable character of this 
complication. 

The Schwinger-Keldysh real-time finite-temperature field-theoretical de- 
scription [jT3|, [T3j] of statistical systems based on the Kubo-Martin- Schwinger 
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(KMS) 07] boundary condition for a field: 



This formal trick introduces into formalism the temperature T = 1/(3 but, 
without fail, leads to the equilibrium fluctuation-dissipation conditions 
IB) (see also [jigfl). Beside this we should have the two-temperature theory 



describing kinetic energy dissipation process (for initial and final states 
separately). It is evident that in such theory with two temperatures it is 
impossible to use the KMS boundary condition. 

In the (S-matrix approach finite-temperature description can be intro- 
duced (e.g. j2Uj and references cited therein) taking into account that, for 
instance, 

dY n = |a n , m |]I /» f 3 oV v e fa) = (9 2 + m /D 1/2 > 
i {Z7c) 6 2e{q t ) 

is the differential measure of final state. Then we can define the temper- 
ature as the function of initial energy through the equation of state, i.e. 
proportional to the mean energy of created particles. Such introduction 
of temperatures as the Lagrange multiplier is obvious for microcanonical 
approach [jnj. The initial-state temperature will be introduced by the 
same way. Using standard terminology [pifl , we will deal with the 'me- 



chanical' perturbations only |2I| and it will not be necessary to divide the 



perturbations on 'thermal' and 'mechanical' ones [^Bj . 

Introducing temperature as the Lagrange multiplier we should assume 
that the temperature fluctuations are small (Gaussian). In opposite case 
the notion of temperature looses its sense. The 'working' idea concerning 
nonequilibrium processes based on assumption that evolution of a system 
goes through few phases. In the first 'fast' phase the s-particle distribu- 
tion functions D s , s > 1, strongly depends from initial conditions. But 
at the end of this phase the system forgets the initial-state information. 
Second phase is the 'kinetic' one. One can expect that the space-time fluc- 
tuations of thermodynamical parameters in this phase are large scale, i.e. 
there is macroscopical domains in which the subsystems are equilibrium, 
with Gaussian fluctuations of thermodynamical parameters. In the last 
'hydrodynamical' phase the whole system is described by macroscopical 
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parameters. We will see that the Schwinger-Keldysh [JT3|, [T4], [TTj formalism 
is applicable for 'hydro dynamical' phase only. 

The above described ^-matrix finite-temperature description can be re- 
alized not only for uniform temperature distribution (we have done first 
step in this direction wishing to introduce initial and final temperatures 
separately). So, introducing cells of measuring device (calorimeter) and 
introducing the energy-momentum shells of each cell separately we can in- 
troduce the individual temperatures in each cell. This can be done since in 
the S'-matrix theory the measurement performed by free (mass-shell) par- 
ticles, i.e. the measurement of energy (and momentum) can be performed 
in each cell separately. This allows to capture the 'kinetic' phase also (if 
the number of calorimeter cells is high enough). In this phase multiparticle 
distribution functions D s , s > 1, are functionals of one-particle distribution 
function Di only. This means the 'shortened' description of the nonequi- 
librium medium J25j. We will return to this question in Sec. 4 considering 
the dissipative processes thermal descriptions applicability range. 

The microcanonical description assumes that the energy of system is 
known with arbitrary accuracy. Introducing the measurement cells and 
corresponding energy shells we assume that the energy in each cell can be 
measured with arbitrary accuracy. It is why we should work in the frame 
of Wigner functions formalism [j25[ . 



E. Wigner had offered the function W(q, R) for the quantum states 
phase space description 



W{q, R)= J dre iqr ^{R + r/2)V*(R - r/2), 

where ^f(x) is the wave function of state. The existence of other approaches 
must be mentioned [P7fl . But as will be seen below the Wigner's description 
is mostly natural for us. 

In the classical limit h = the function W(q, r) coincides with the phase 
space probability distribution function. It obeys the equation |28j]: 



W = {W,H} + 0(h), 

which coincides with Liouville equation only in the classical limit % = 0. 

The extension of Wigner's idea on the relativistic case uses the con- 
nection between Wigner's approach and inclusive description of inelas- 
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tic scattering processes [|25|, pgfl . But the Wigner functions are not di- 
rectly measurable quantities because of the quantum uncertainty principle 
AqAr ~ U. Just this restriction leads to impossibility to take the measure- 
ment (calorimeter) cells 4-dimension Ar arbitrary small and defines the 
natural boundary of Wigner functions approach applicability. Wishing to 
use the Wigner-functions description of experiments the corresponding the- 
ory must take into account this restriction. The discussion of this question 
is given in Sec. 4. 

So, in our terms one can use the thermodynamical formalism if the non- 
stationary mediums 'shortened' description may be applied: in this case 
a mean value of correlation functions over the space-time are negligible 
and the fluctuations of thermodynamical parameters are small (Gaussian). 
Other approach should be mentioned also. Proposing |32j that the equilib- 
rium in the nonuniform nonstationary medium may be attained in small 
regions more quickly than in the whole system the entropy maximalness 
in this restricted domains of a system can be used for construction of the 
'local equilibrium density matrix' (LDM) But LDM is applicable 

for description of processes in which dissipation may be disregarded [f3Uj. 
Nevertheless, if the energy-momentum density of nonstationary flow is con- 
siderably smaller than the energy density of matter then first one can be 
taken into account perturbatively considering LDM as the initial condition. 
This modifies LDM to 'nonstationary density matrix' (NDM) of Zubarev 
P2| introducing an infinitesimal interaction with a heat bath to get the 
increasing entropy. We will return to this question in Sec. 5. 

The /S-matrix will be introduced phenomenologically, using ordinary in a 
quantum field theory reduction formalism. This leads naturally to necessity 
to introduce the boundary conditions for interacting fields $((100), where 
doc is the infinitely far hypersurface, e.g. The value of ^((Joq) specify 
the environment of a system. 

We start from vacuum boundary condition $(croo) = familiar for a 
field theory. This theory can be applied in the particle physics. The sim- 
plest choice of $(cToo) 7^ assumes that the system under consideration 
is surrounded by black-body radiation. Just this 'boundary condition' re- 
stores the Schwinger-Keldysh [jTTj real-time finite-temperature field theory 
T2f from ^-matrix formalism in the 'hydrodynamical' phase and gives the 
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dynamical interpretation of the KMS periodic boundary condition. 

One should admit also that last choice of boundary condition is not 
unique: one can consider another organization of the environment of con- 
sidered system. The ^-matrix interpretation is able to show the way of 
adoption of formalism to the arbitrary environment^]. It should broaden 
the potentialities of the real-time finite-temperature field-theoretical meth- 
ods. For instance, for heavy nucleus high energy interactions. The special 
interest represent also the topological effects, but, by above mention reason, 
in this paper consideration will be performed in the perturbation theory 
framework only (see also concluding Section). 

The central purpose of this review paper is to describe connections be- 
tween ordinary ^-matrix description and popular in the modern literature 
real-time finite-temperature field theories. We wish to discuss: 
-The QCD jets dominance in deep asymptotics over n (Sec. 2); 
In this section we would like to show why the real — time formalism is 
needed for our dissipative process description. 

-The ^-matrix interpretation of Schwinger-Keldysh theory (Sec. 3); 
In this section the uniform temperature description of the state will be 
introduced in the spirit of microcanonical description. It is shown the way 
of explicit calculations to show coincidence of used microcanonical descrip- 
tion and ordinary (Gibbs) canonical formalism. 

-The range of applicability of finite-temperature description (Sec. 4); 

In this section the necessity and sufficiency of Bogolyubov's 'shortened' 

description is discussed. 

-The ^-matrix description of media with nonuniform temperature distri- 
bution (Sec. 5); 

In this section the Wigner functions formalism is introduced. The range if 
its applicability to describe an experiment is shown. 

-The comparison of our ^-matrix approach with 'nonstationary statistical 
operator' of Zubarev |32j (Sec. 6). 

In this section the main distinction between ^-matrix (microcanonical) and 
Zubarev's (canonical) perturbation theories is shown. 
-Concluding remarks (Sec. 7). 

3 This question was considered also in |29| 
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In this section the way as the nonperturbative effects may be included in 
the formalism is discussed. 



2 Phenomenology 



To build the phenomenology |33j of high multiplicity processes let us intro- 



duce the classification of asymptotics over n. For this purpose it is useful 
to consider the 'big partition function': 

T ( z , s ) = J2 z n a n {s), T(l, s) = atot(s). 

n 

Strictly speaking, summation over n is performed up to n max . But we 
can extend summation up to infinity^ if the weight z is sufficiently small, 
< z < z max . So, T(z, s) can be considered as the nontrivial function of z 
with sufficient accuracy. Note that z max > 1 since a n (s) decrease with n. 

If we know T(z : s) then cr n (s) is defined by inverse Mellin transformation. 
This gives (usual in thermodynamics) equation (of state): 

d 

n = z—lnT{z,s) (2.1) 

Solving this equation we can estimate the asymptotics of o n \ 

a n (s) ~ e - nln ~ zM , (2.2) 



where 1 < z(n,s) « z max is smallest solution of eq. QZTT 



It follows from fl2Tg ) that at n — > oo the solution of (|ZT| ) must tend to 



singularity z s of T(z, s) and the character of singularity is not important. 
So, we must consider three possibility: 

a). z s = z a = 1, b). z s = z b = oo, c). z s = z c , 1 < z c < oo. 



Following to Lee and Yang [33j there is not singularities at < z < 1. 

Let us consider now the physical content of this classification, 
a) z s = 1. 

It is known that the singularity z s = 1 reflects the first order phase transi 
tion M. To find a n for this case we would adopt Langer's analyses 



l I.e. wishing to consider the 'thermodynamical limit'. 
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Introducing the temperature 1/(3 instead of total energy yfs we can use the 
isomorphism with Ising model. For this purpose we divide the space vol- 
ume on cells and if there is particle in the cell we will write (-1). In opposite 
case (+1). It is the model of lattice gas well described by Ising model. We 
can regulate the number of down-looking spins, i.e. the number of created 
particles, by the external magnetic field H. Therefore, z = exp{— /3H} and 
H is the chemical potential. 

The corresponding partition function in the continuous limit (see 
also [001) has the form: 



R(P,z) = J Dfie-^ dx ^ 2 - £ ^ +afl4 - x ^, (2.3) 

where e ~ (1 — and A ~ H, with critical temperature 1/(3 C . 

If j3 c > (3 there is not phase transition and the potential has one min- 
imum at ji = 0. But if j3 c < j3 there is two degenerate minimum at 



fi± = ±Je/2a if A = 0. Switching on H < the left minimum at 



^_ ~ —ye/2a becomes absolute and the system will tunnel into this mini- 
mum (see also |57|| ). This process describes particles creations as a process 
of spins trippings. 

The eq . flZTT] ) gives at n — > oo 

\nz ~ n~ 1/3 > 0. 

In result, 

a n ~e- an2/3 >0(e- ?1 ), a>0, 

i.e. decrease slower then e~ n . The quasiclassical calculation shows that the 
functional determinant is singular at H = 0. It must be underlined that 
in the used Ising model description the chemical potential deforms ground 
state. In result the quasiclassical approximation is applicable since In z << 
1, i.e. since the processes of spin trippings are rear at high multiplicity 
region. It is easy to show in this approximation [j315j that the functional 
determinant is singular at H = 0, i.e. at z = 1. 

Note that z decrease to one with n. This unusual phenomena must 
be explained. Considered above mechanism of particles creation describes 
'fate of false vacuum' [|37j| . In the process of decay of unstable state the 



clusters of new phase of size X are created. If the cluster have dimension 
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X > X c its size increase since the volume energy (~ X 3 ) of the cluster 
becomes better then the surface tension energy (~ X 2 ). This condition 
defines the value of X c . The 'critical' clusters wall accelerate, i.e the work 
needed to add one particle into cluster decrease with X > X c . This explains 
the reason why z decrease with n noting that In z is proportional to Gibbs 
free energy per one particle. 

The described mechanism of particles creation assumes that we had 
prepared the equilibrium system in the unstable phase at ji + ~ + Je/2a 



and going to another state at /i_ ~ —ys/2a the system creates particles. 
The initial state may be the QGP and final state may be the hadrons 
system. Therefore, we must describe the way as the quarks system was 
prepared. 

Following to Lee- Yang's picture of first order phase transition |3| 



see 



also there is not phase transition in a finite system (the partition func- 



tion can not be singular for finite n max ). This means that the multiplicity 
(and the energy) must be high enough to see described phenomena, 
b) z s = oo. 

Let us return to the integral (|2T3f) to investigate the case (3 C > (3. In this 



case the potential has one minimum at /i = 0. The external field H cre- 
ates the mean field Jl = /i(H) and the integral (^3) should be calculated 
expanding it near fi = ft. In result, in the quasiclassical approximation (ji 
increase with increasing n), 

]nR(J3,z) ~ (lnz) 4/3 . 

This gives Inz ~ n 3 and In a n ~ — n 4 , i.e. 

on < 0(e~ n ). 

There is also other possibility to interpret considered case b). For this 
case we can put 

lnT(z, s) = n (s) + n(s)(z - 1) + 0((z - l) 2 ) (2.4) 

at \z—l | << 1. By definition no (s) = lna t ot- The experimental distribution 
of In T(z,s) — no(s) for various energies shows that the contributions of 
0((z— l) 2 ) terms increase with energy j38j. The hadrons 'standard model' 
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(SM) assumes that 



\nt(z, s) = uq(s) + n(s)(z — 1) 

is the Born term in the perturbation series ( |2.4|) . There is various interpre- 
tations of this series, e.g. the multiperipheral model, the Regge pole model, 
the heavy color strings model, the QCD multiperipheral models, etc. In 
all this models no = a\ + a<i In s, < a<i « 1 and n(s) = b\ + 62 Ins, 
62 > 0. The second ingredient of hadrons SM is the assumption that mean 
value of created particles transfers momentum < k >= const, i.e. is the 
energy (and multiplicity) independent. It can be shown that under this 
assumptions: 

In T(z, s) = n {s) + £ c n {s){z - If, c x = n (2.5) 

n 

is regular at finite values of z and is able to give well confirmed by 



experiment predictions. 

Inserting (|2T51 ) into ( |2TT|) we find that z(n, s) is the increasing function 
of n. Therefore, 

o n < 0(e~ n ). (2.6) 

But the SM have a finite range of validity: beyond n ~ n 2 the model 
must be changed since it is impossible to conserve < k >= const, at higher 
multiplicities . 



We should underline once more that only two possibilities a) and b) can 
be deduced from representation ( PT3|) , see also [P5fl . But nevertheless there 
is other possibility: 
c) 1 < z s < 00. 
Let us assume now that 

T(z, s) ~ (1 - ^rY\ 7>0. (2.7) 
z c 1 

Then, using normalization condition, (dT(z, s)/dz)\ z= \ = n\j{s) we can 
find that z c (s) = 1 + ^f/flj(s). The singular structure ( |2.Y|) is impossible 



in SM because of condition < k >= const. But if \z — 1| << 1 we 
have estimation ( \Z.4Jj ). The difference between SM and c) is seen only at 



1 — (z — l)/{z c — 1) << 1, i.e. or in asymptotics over n or in asymptotics 
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over energy. The singular structure is familiar for 'logistic' equations of 
QCD jets, e.g. [HQ . 



In considered case z = z c + 0(nj/n) and at high energies (rij(s) >> 1) 

a n ~ e" 7n/ ^ = 0(e~ n ). (2.8) 

Therefore, 

comparing ( \2.tJj) and j\2. we can conclude that at sufficiently high ener- 



gies, i.e. if n~j >> n, where n is the SM mean multiplicity, the mechanism 
c) must dominate in asymptotics over n. 

It is the general, practically model independent, prediction. It has im- 
portant from experimental point of view consequence that at high energies 
there is wide range of multiplicities where the SM mechanism of hadrons 
creation is negligible. In other words, the CQGP of high multiplicity pro- 
cesses is the dynamical consequence of jets and SM mechanisms. At tran- 
sition region between 'soft' of SM and 'hard' of jets one can expect the 
'semihard' processes of minijets dominance. 

The multiplicity distribution in jets has interesting property noted many 
decades ago by Volterra in his mathematical theory of populations [BTJ • In 



our terms, if one-jet partition function has the singularity at z]p{s) = 
1 + r y/nj(s) then two-jet partition function must be singular at 

and so on. Therefore, at high energies and n > Uj(s) the jets number must 
be minimal (with exponential accuracy). This means that at n — > oo the 
processes of hadrons creation have a tendency to be Markovian (with sharp 
increase of transverse momentum < k >) and only in the last stage the 
(first order) phase transition (colored plasma) — ► (hadrons) may be seen. 

One can say that in asymptotics over n we consider the 'inflational' 
chanel of thermalization which is so fastQ that the usual confinement forces 
are 'freezed' and do not play important role in final colored plasma creation. 

5 The partons life time with virtuality \q\ is ~ l/\q\ and the time needed for hadrons of mass rrih 
formation is ~ 1/tto/,. Therefore the parton have a time to decay before hadrons formation if \q\ » m. 
But this situation is rear since the thermal motion in the initial stage of process is high. 
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3 5-matrix interpretation 



3.1. Vacuum boundary conditions. 

The starting point of our calculations is n- into m-particles transition 
amplitude a„ jm , the derivation of which is well known procedure in the 
perturbation theory framework. For this purpose the (n + m)-point Green 
function G n m is introduced fS2|l . To calculate the nontrivial elements of 



S'-matrix one must put the external particles on the mass shell. Formally 
this procedure means amputation of the external legs of G c n m and further 
multiplication on the free particles wave functions. In result the amplitude 
of m- into n-particles transition a nm in the momentum representation has 
the form: 

fln,m((?)»l (P)m) = H)^™ II 0(<fc) II & fa)Z (</>) . (3.1) 

k=l k=l 

Here we introduce the 'annihilation' operator 

0(g) = / dxe-^(x), fa) = j^- y (3.2) 

(jf(pk) is the 'creation' operator and qt and pk are the momentum of in- 
and out-going particles. In (|3.1| ) 



Z(<f>) = J D$e lS{<p) - 



■iV(G>+<j>) 



is the generating functional. The total action was divided on two parts, 
where 5(3>) is the free part and V(<&, 4>) describes the interactions. At the 
very end one should put the auxiliary field = 0. 

To provide the convergence of the integral ( |3.1|) over scalar field $ the 
action must contain positive imaginary part. Usually for this purpose 
Feynman's ie-prescription is used. It is better for us to shift infmitesimaHy 
time contour to the upper half plane [|3], [IUj], i.e. to the Mills contour 



C + : t — > t + ie, e > 

and after all calculations to return the time contour on the real axis, e 
+0. 
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In eq. (|3.I|) the integration is performed over all field configurations 
with standard vacuum boundary condition: 



/ (fxd^d^) = I da^d^ = 0, 



which assumes zero contribution from the surface term. 

Supposing that the particles number and momenta are insufficient for 
us we introduce the probability 



|2 



1 r n n 

r(P) = E -r-T / du n {q)du m {p)5^\P - £ q k )^\P - £ ft)k, m| ■ 

(3.3) 

where 

n n d^qk 

du n (q) = J] dw{q k ) = JJ (27r) 3 2e(g fc ) ' 6 = ^ + m h) i 

is the Lorentz- invariant phase space element. We assume that the energy- 
momentum conservation ^-function was extracted from the amplitude. It 
was divided onto two parts: 

£ (4) (£^ - Eft) = / d 4 P5^(P -T,q k )6 {4) (P - Eft) (3-4) 
It is not too hard to see that, up to phase space volume, 

r = J d 4 Pr(P) 

is the imaginary part of amplitude < vac\vac >. Therefore, computing 
r(P) the standard renormalization procedure can be applied and the new 
divergences will not arise in our formalism. 

The Fourier transformation of <5-functions in ( |5.3|) allows to write r(P) 
in the form: 



where 



1 71 7JI 

p(«i,« 2 ) = E — j — r / U{dco(q k )e-^} R {cMftK^lKJ 2 . (3.5) 
n ,m n\m\ J k= i k=i 

Introduction of the 'Fourier-transformed' probability p(a\, a?) means only 
that the phase-space volume is not fixed exactly, i.e. it is proposed that 
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4-vector P is fixed with some accuracy if on are fixed. The energy and 
momentum in our approach are still locally conserved quantities since the 
amplitude a nm is translational invariant. So, we can perform the transfor- 
mation: 

Oil Qk = (CKI - (Ti) Qk + Ox Qk -> (Oil - CF\) Qk + 0\P 

since 4-momenta are conserved. The choice of o\ fixes the reference frame. 
This degree of freedom of the theory was considered in J33], [I5[ . 
Inserting (p71|) into ( gT5Q we find that 



Z 1 s\ s\ 

p(ai,oi2) = exp{i / dxdx'((f) + (x)D + ^(x — x ,a2)<fi-(x) — 

-}-(x)D-+(x - x', a 1 )0 + (x , ))}Z(0 + )Z*(0_), (3.6) 

where and D + are the positive and negative frequency correlation 

functions: 

D + 4x - x', a) = -i J duj(q)e iq{x " x '~ a) 

describes the process of particles creation at the time moment xq and its 
absorption at x' Q , xq > x' Q , and a is the center of mass () 4-coordinate. 
Function 

D_+0 -x',a)=ij dco{q)e^ x ~ x ' +a) 

describes the opposite process, xq < x' . These functions obey the homo- 
geneous equations: 

(d 2 + m 2 ) x G + - = (d 2 + m 2 ) x G- + = 

since the propagation of mass-shell particles is described. 

We suppose that Z((p) may be computed perturbatively. For this pur- 
pose following transformation will be used: 

e -iV{4>) = e -i J dxj(x)4>'(x) e i J dxj{x)cj>{x) e -iV \4>') = 

= e Jdx<Kx)&(x) e -%VW) = 

= e -iV(-ij) e ifdxj(x)(t>(x)^ (3-7) 

where (j) was defined in (pT^). At the end of calculations the auxiliary 
variables j, (j)' should be taken equal to zero. Using the first equality in 
(|377p we find that 

_ e -iJdxj(x)$(x) e -iV{<5>+(j)) e -i I Jdxdx'j(x)D ++ (x-x')j(x')^ p^gj 
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where D ++ is the causal Green function: 

{d 2 + m 2 ) x G ++ {x - y) = 5(x - y) 

Inserting ( |5.8|) into (|5.fci| ) after simple manipulations with differential oper- 
ators, see ( |3.Y| ) we find the expression: 

p(a h a 2 ) = e~ iV ^ +iV ^ x 
x exp{- / dxdx'(j + (x)D + -(x — x\ ai)j-(x') — 

j-{x)D_ + (x - x' : a 2 )j+(x') - 
-j + (x)D ++ (x - x')j+(x') + j.(x)D—(x - x')j-(x'))}, (3.9) 

where 

L>__ = (D ++ Y 

is the anticausal Green function. 

Considering the system with large number of particles we can simplify 
calculations choosing the CM frame P = (Pq = E,0). It is useful also 
20[ , to rotate the contours of integration over ao^'- «o,fc = — i>Pk, Imfik = 



0, k = 1,2. In result, omitting unnecessary constant, we will consider 
p = p((3 u (3 2 ). 

External particles play the double role in the ^-matrix approach: their 
interactions create and annihilate the interacting fields system and, on 
the other hand, they are probes through which the measurement of the 
system is performed. Since (3k are the conjugate to the particles energies 
quantities we will interpret them as the inverse temperatures in the initial 
(/?+) and final (/?_) states of interacting fields. But there is the question: 
are constants /?& really the 'good' parameters to describe the system. 

The integrals over /?&: 

r(E) = [ ^*/k e & + M E e- F &>M, (3.10) 

where 

F(p 1 ,[3 2 ) = -lnp((3 1 ,(3 2 ), 

can be computed by the stationary phase method. This assumes that the 
total energy E is a fixed quantity. The solutions of the equations (of state): 

E-°*$M, k-1,2, (3.1!) 
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gives the mostly probable values of (5k at a given E. Eqs. always 
have the real solutions and, because of energy conservation law, both eqs. 
(jXTTj) have the same solution with the property [jltjj] : 



(5 k = (5(E), (5>0. 

Assuming that (5 is the 'good' parameter, i.e. the fluctuations of (5k are 
Gaussian we can interpret F((5i,(52) as the free energy and l/(3k as the 
temperatures. Such definition of thermodynamical parameters is in a spirit 
of microcanonical description. We will return to this question in Sec. 4. 

The structure of generating functional (^) is the same as the generat- 
ing functional of Niemi-Semenoff |T3| have. The difference is only in the 



definition of Green functions which follows from the choice of boundary 
condition (|2T6| ) . The Green functions Dij,i,j = +, — were defined on the 
time contours C± in the complex time plane (C_ = C+). This definition of 
the time contours coincide with Keldysh' time contour [M. The expres- 



sion can be written in the compact form if the matrix notations are 
used. Note also a doubling of the degrees of freedom. This doubling is 
unavoidable since Green functions Dij are singular on the light cone. 
3.2. Closed-path boundary conditions. 

The generating functional p(ai,ot2) has important factorized structure, 
see (W^i)' 

p(a 1 ,a 2 ) = e^ a ^p (<f> ± ), 

where the operator 

N(ai, a% 4>) = J dxdx'((j) + (x)D + -(x — x', a2)4>~(x') — 

-(j)-{x)D_ + (x - x\ ai)j) + (x')) (3.12) 

acts on the generating functional 

Po{<f> ± ) = Z{<f> + )Z*((f>^ = 
= J D^ + D^^e lS ^ + ^ tS ^-^ iV ^ ++ct)+)+tV ^- +(f) -\ (3.13) 

of measurables. All 'thermodynamical' information was contained in the 
operator N(ai,ot2', 4>) and interactions are hidden in po(4>±). One can say 
that action of the operator N maps the system of interacting fields on the 
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measurable states. Last ones are 'labeled' by a\ and oli- Just this property 
allows to say that we are dealing with 'mechanical' fluctuations only. To 
regulate the particles number we can introduce into N dependence from 
'activities' Z\ and z<i for initial and final states separately. 

The independent fields 0+, 0_ and $+, were defined on the time con- 
tours C + , C-. By definition, path integral ( [3.130 describes the closed path 
motion in the space of fields We want to use this fact and introduce a 
more general boundary condition which also guaranties cancelation of sur- 
face terms in the perturbation framework. We will introduce the equality: 

/ da^+d^+ = I dafl-d"®-. (3.14) 

The solution of eq. (|37T^) requires that the fields $ + and <3>_ (and theirs 
first derivatives d^±) coincide on the boundary hypersurface a^: 

$±(^oo) = $(<7oo), 

where, by definition, <I>(<7oo) is the arbitrary, 'turning-point', field. 

The existence of nontrivial field $(croo), in absence of surface terms, has 
influence only on the structure of Green functions 

G++ =< T<$> + <$> + >, G+ =< >, 

G + =<$_$+>, G =<f $_$_>, (3.15) 

where T is the antitemporal time ordering operator. This Green functions 
must obey the equations: 

(d 2 + m 2 ) x G + 4x — y) = (d 2 + m 2 ) x G„ + (x - y) = 0, 
(d 2 + m 2 ) x G ++ (x -y) = (d 2 + m 2 ) x G-4x - y) = 5(x - y), (3.16) 

and the general solution of these equations: 

Ga Dn -\- ya-, 
Gij = gij, i^j (3.17) 

contain the undefined terms which must obey the homogenous equa- 
tions: 

(d 2 + m 2 ) x9ij (x-y) = 0, U = +,-. (3.18) 
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The general solution of these equations (they are distinguished by the 
choice of the time contours C±) 

9ij(x - x') = / dLu(q)e^ x -^n l3 (q) (3.19) 

are defined by the functions riij. Last ones are the functionals of 'turning- 
point' field ^((Joo): if $(coo) = we must have = and we will come 
back to the theory of previous section. 

Our aim is to define n^. We can suppose that 

n tj ~< $(<7oo) • • • $(<7oo) > • 

The simplest supposition gives: 

riij ~< Q&j >~< $ 2 (<t°°) > . (3.20) 

We will find the exact definition of starting from the ^-matrix inter- 
pretation of the theory. 

We should suppose there are only free, mass-shell, particles that on the 
infinitely far hypersurface doo. Formally this follows from (|5.lY|) -( gTT9| ) 



and is natural in the ^-matrix framework [|7[. In other respects choice of 
the boundary condition is arbitrary. 

Therefore, our aim is the description of evolution of the system in a 
background field of mass-shell particles. We will assume that there are 
not any special correlations among background particles and will take into 
account only the energy-momentum conservation laws constraints. Quan- 
titatively this means that multiplicity distribution of background particles 
is Poison-like, i.e. is determined by the mean multiplicity only. This is in 
spirit of definition of in eqs.( |3.19|) , (|3.2U|) . 



Our derivation is the same as in |55[]. Here we restrict ourselves men- 



tioning only the main quantitative points. 

In the vacuum case of Sec. 3.1 the process of particles creation and theirs 
further absorption was described. In presence of the background particles 
this time-ordered picture is wiped out: appears possibility of particles ab- 
sorption before theirs creation. 

The particles creation and absorption was described by the product of 
operator exponent (|3T6|) . One can derive (see also |fo5(| ) the generaliza- 



tions of (|53): presence of the background particles will lead to the same 
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structure: 



Pcp = e «(«*,) Po (^ : 



where Rq(4>±) is the same generating functional, see (|3.13|) . But the oper- 
ator N((/)*(j)j), i, j = +, — , should be changed wanting to take into account 
the external particles environment. 

The operator (j)*(q) was interpreted as the creation and <pi(q) as the 
annihilation operator, see definition (|3.1|). Correspondingly the product 
4>*{q)4>j{q) acts as the activity operator. So, in the expansion of N{<j>*<f)j) 
we can leave only first nontrivial term: 

NWh) = J du(q)i>:(q)n tJ ^(q), (3.21) 

since no special correlation among background particles should be ex- 
pected. If the external (nondynamical) correlations are present then the 
higher powers of will appear in expansion ( [3.21|) [EZ9f . Following to 



the interpretation of we conclude that is the mean multiplicity of 
background particles. In (|3.21| ) the normalization condition N(0) = was 



used and summation over all i,j was assumed. (In the vacuum case only 
the combinations i ^ j was present.) 

Computing R cp we must conserve the translational invariance of am- 
plitudes and extract the energy-momentum conservation (5-functions. We 
must adjust to each vertex of in- going particle in <2 n>m the factor e - JQ W 2 
and for each out-going particle e~ ia2q ^ 2 one after Fourier transformation of 
this (5-functions. 

So, the product e -* Q W 2 e -* a j<?/ 2 can De interpreted as the probability 
factor of the one-particle (creation^ annihilation) process. The n-particles 
(creation^ annihilation) process' probability is the simple product of these 
factors if there is not special correlations among background particles. This 
interpretation is evident in the CM frame ct^ = (— 0). 

After this preliminaries it is not hard to find that in the CM frame we 
have: 

n ++ {q ) = n__(g ) == = fi^ f 1 + ^ 2 ). (3.22) 

Computing n^ for i ^ j we must take into account that we have one 
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additional particle: 

n+-(do) == 6(g )(l + ftWi)) + B(-g )n(-g A) (3.23) 

and 

ra-+(«b) = B(go)nM 2 ) + 6(-gb)(l + n{-q (3 2 )). (3.24) 

Using (\6.2'4) , (|3.^3|) and ( |5.^4|) , and the definition ( |3.lY| ) we find the Green 
functions (the matrix Green functions in the real-time finite-temperature 
field theories was introduced firstly in 

G«(x - x', 09)) = / (0i^<"'>G j3 -( g , 09)) 



where 







q^—m—te 



1 n(Ako|)a-(A) »(^l5b|) 1 



and 

a±(/?) = -et< 

The corresponding generating functional has the standard form: 

Pc P {j±) = exp{-iV(-ij+) + iV(-ij_)} x 
x exp{^ J dxdx' ji(x)Gij(x - x', ((3))jj(x')} (3.26) 

where the summation over repeated indexes is assumed. 

Inserting ( |3.26|) in the equation of state ( |3.11| ) we can find that (3\ = 
= (3(E). If (3(E) is a 'good' parameter then Gij(x — x'; (3) coincide with 

the Green functions of the real-time finite-temperature field theory and the 

KMS boundary condition: 

G + _(t - t') = G-+(t - 1' - i/3), G-+(t - t') = G+-(t -t' + i/3), (3.27) 

is restored. The eq. (|3727| ) can be deduced from (pT25|) by the direct calcu- 
lations. 
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4 Applicability of finite-temperature description 



4-1. The Schwinger-Keldysh formalism. 

There is various approaches to build the real-time finite-temperature 
field theories of Schwinger-Keldysh type (e.g. iPDjl ). All of them uses 



various tricks for analitical continuation of imaginary-time Matsubara for- 
malism to the real time p7j| . The basis of the approaches is introduction 
of Matsubara field operator 

$ M (x,/3) = e^$ 5 (xe-/^, (4.1) 

where ^(x) is the interaction-picture operator, instead of Heisenberg op- 
erator 

$(x,£) =j tH $s(xe- itH . 

This introduces the averaging over Gibbs ensemble instead of averaging 
over zero-temperature vacuum states. 

If the interaction switched on adiabatically at the instant ti and switched 
off at tf then there is the unitary transformation: 

$(x) = U(ti,t f )U(U,t)$s(x)U(t,ti)' 

Introducing the complex Mills time contours |E3j] to connect U to t, t to tf 
and tf to ti we form 'closed-time' contour C (the end-points of the contours 
C + and C_ are joint together). This allows to write last equality in the 
compact form: 

where Tq is the time-ordering on the contour C operator. 

The corresponding expression for the generating functional Z(j) of cor- 
relation (Green) functions has the form: 

Z(j^) = i?(0) < Tce i fc d4x { L ir>.t{x)+j(x)^(x)}s > 

where <> means averaging over initial state. 

If the initial correlations have little effect we can perform averaging over 
Gibbs ensemble. This is the main assumption of formalism: the generating 
functional of the Green functions Z(j) has the form in this case: 



Z(j) = [D& < &;t i \e-P H T c e i fc d ^( x M x )\& ; t i 



> 
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with $' = $'(x). In accordance with Q4.1|) we have: 



and, in result, 

Z(j) = J D&e !c ? ^W+^W*) (4.2) 
where path integration performed with KMS periodic boundary condition: 

<$>{U) = <$>{U-i(3). 



In (fO|) the contour Cp connects t{ to tf, tf to t{ and to U — if3. Therefore 



it contains imaginary-time Matsubara part ti to ti — i(3. More symmetrical 
formulation uses following realization: ti tot/, tf to tf — i/3/2, tf — i(3/2 
to ti — i/3/2 and t{ — i(3/2 to t{ — i(3 (e.g. |T2j]). This case also contains the 



imaginary-time parts of time contour. Therefore, eq. ( |4.4 ) presents the 
analitical continuation of Matsubara generating functional to real times. 

One can note that if this analitical continuation is possible in Z(j) 
then representation (|Oj) gives good recipe of regularization of frequency 
integrals in the Matsubara perturbation theory, e.g. jKj], but nothing 



new for our problem since the Matsubara formalism is a formalism for 
equilibrium states only. 

Taking ti = — oc and tf = +00 and calculating integral ( |4.2|) perturba- 



tively we find coincidence of Z(j) and R(/3) from (|3.2(d| ) with Green func- 



tions defined in ( |3.25|) if f3\ = fy. This 'factorization' of contributions from 



contours C + and C_ in the integral ( |4.2|) follows from Rieman-Lebesque 
lemma [^] which is applicable in the perturbation framework [^3| , [T2J] . 



Note absence of Matsubara parts of contour, which prevents the factoriza- 
tion, in the derived '^-matrix generating functional' ^tj|) by definition 
(importance of this circumstances is discussed in Sec. 7). 

4-2. Range of the 'hydrodynamical 7 approximation. 

Let us return now to eq.( |5.h|) . To find the physical meaning flip) we must 
show the way as they can be measured. If there is nonequilibrium flow it 
is hard to invent a thermometer (or thermodynamical calorimeter) which 
measures locally in space-time the temperatures of this dissipative pro- 
cesses. But there was described another way - to define the temperatures 
through equations of state. This is possible in the accelerator experiments 
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where the total energy E is fixed. So, we will define fii^) through equations 
of state i.e. considering 1/A(2) as the mean energy of particles in 

the initial (final) state. But even knowing solutions of this equations one 
can not find p(E, z) correctly if the assumption (3u2) are 'good' quantities 
is not added, i.e. that the fluctuations near solutions of eqs.( |37TT|) are small 
(Gaussian). 

This assumption is the main problem toward nonequilibrium thermo- 
dynamics. The problem in our terms looks as follows: the expansion near 
(3\{2){E) gives asymptotic series over 

/ D .s ~ / U{du(k t )dr l } < e(k 1 )e(k 2 ) ■ ■ ■ > |( ri ,r 2 ,...), 

where <>q means averaging over fields drown on fixed points of phase 
space (/c, r)i. In other words, the fluctuations near [3i(2)(E) are defined by 
value of inclusive spectra familiar in particles physics. Therefore, f3u2)(E) 
are 'good' quantities if this inclusive spectra are small. But this is too 
strong assumption. More careful analysis shows that it is enough to have 
the factorization properties [£E| : 

/ \{{duj{k i )dr i } < e{ki)e{k 2 ) •••> |( n ,r 2 ,...) ~ 

II / dn(ki)dn < e{ki) > | (Pi) ~ 0. (4.3) 

It must be noted that this is the unique solution of problem since the 
expansion near ^u 2 ){E) unavoidably leads to asymptotic series with zero 
radii of convergence. 

One can hope to avoid this problem working permanently in the energy- 
momentum representation, i.e. without introduction of temperatures. Of 
course this is possible in particles physics, but if (3u2)(E) is not the 'good' 
parameter this means that all correlations between created particles are 
sufficient, i.e. only the energy-momentum representation did not solve the 
problem. 

At the end, discussed factorization property of D s , s > 1, is well known 
Bogolyubov's condition of 'shortened' description of nonequilibrium ther- 
modynamical systems with s-particle distribution functions D s , s > 1, 
expressed in terms of Di. It is the condition for the 'hydrodynamical' de- 
scriptions applicability since it assumes that the constant (3\{E) = (3 2 {E) 
is a 'good' parameter for description of whole system. 
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Considering a problem with nonzero nonequilibrium flaw it is hard to 
expect that j3u2\{E) is a good parameter, i.e. that the factorization condi- 
tions are hold. Nevertheless, as was mentioned above, there is possibility 
to have the mean values of correlators sufficiently small in restricted ranges 
of phase space. It is the so called 'kinetic' phase of the process when the 
memory of initial state was disappeared, the 'fast' fluctuations was aver- 
aged over and we can consider the long-range fluctuations only. 



5 Local equilibrium hypothesis 

Let us return now to description of experimental situation in the high 
multiplicity experiments. Having at energies of modern accelerators thou- 
sands of particles in a final state it is a hard problem even to count such 
big numbers. So, the number of particles n can not be considered as a 
trigger. Moreover, it seems naturally that it is not important have we hun- 
dred thousand of particles or hundred thousand plus one. To do first step 
toward CQGP it is enough to be sure that on experiment the transition 
of 'hot' initial state into 'cold' final one is examined. For this purpose the 
ordinary calorimeters can be used j5T| . 



So, we must assume that the energies of created particles E{ < £o, where 
£o is fixed by experiment. Then using energy conservation law at given £q 
the number of created particles is bounded from below: n > y/s/eo = n m i n . 
With this constraint the integral cross section 

is measured. Choosing n m i n >> h, i.e. £o << Vs/h(s), we get into high 
multiplicity region. There is also a theoretical possibility to restore the 
quantity ~ o n calculating the difference cr £o (s) — a £o+ se { s ) IMD • 



It is not necessary to measure energy of each particle to have n m i n >> h. 
Indeed, let £{ is the energy of z-th group of particles, e\ + £2 + ••• + £k = ^/s 
and let hi is the number of particles in the group, hi + h^ + ... + hk = ^Q- 
Then, if §i < £q, i = 1, 2, k, we have inequality: k > n min . Therefor, we 



6 It is assumed that the number of calorimeter cells K > k. 
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get into high multiplicities domain since n > k, if £q << i/s/n(s). We can 
use the calorimeter demanding that the educed in each cell energy e\ < £o. 

The preparation of such experiment is not hopeless task and it may 
be sufficiently informative. This formulation of experiment we will put 
in basis of the theory. Theoretically we should shrink the 4-dimension of 
calorimeter cells up to zero since we do not know ad hoc the cells dimension. 
Then the cells index i is transformed into the position of particle r. So we 
come to contradiction with quantum uncertainty principle. This forces to 
use the Wigner functions formalism and the first question which must be 
solved is to find a way as this formalism can be adopted for description of 
our experiment (there is also interesting ideas concerning applicability of 
Wigner functions in fl5TH ). 

5.1. Vacuum boundary condition. 

We start consideration from the assumption that the temperature fluc- 
tuations are large scale. In a cell the dimension of which is much smaller 
then the fluctuation scale of temperature we can assume that the tem- 
perature is a 'good' parameter. (The 'good' parameter means that the 
corresponding fluctuations are Gaussian.) 

Let us surround the interaction region, i.e. the system under consider- 
ation, by N cells with known space-time position and let us propose that 
we can measure the energy and momentum of groups of in- and out-going 
particles in each cell. The 4-dimension of cells can not be arbitrary small 
in this case because of the quantum uncertainty principle. 

To describe this situation we decompose (5-functions in ( |3.4| ) on the 
product of (N + 1) ^-functions: 

n N n v N 

^(P - E Qk) = / II {dQAQu - E q k ,,)}S^(P - E Qu), 

k=l v=\ k=l v=\ 

where are the momentum of k-th in-going particle in the z^-th cell 
and Q v is the total 4-momenta of n v in-going particles in this cell, v = 
1, 2, N. The same decomposition will be used for the second ^-function 
in (|5.4| ). We must take into account the multinomial character of particles 
decomposition on TV groups. This will give the coefficient: 

n\ N ml N 

-5 K {n-Y. n v) — ; — ' r^(m-^m^), 



n x \---n N \ v= i m 1 \---m N \ v=x 
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where 5k is the Kronecker's 5- function. 
In result, the quantity 



r((QU(P)N)= E / \Hn, 



l 2 x 



(n.m) 



x n { n ^^\q„ - 1 n - "± »„» ^ 

u=l k=l n v- k=l k=l m v- k=l 

describes a probability to measure in the z^-th cell the fluxes of in-going 
particles with total 4-momentum Q v and of out-going particles with the 
total 4-momentum P v . The sequence of this two measurements is not fixed. 
The Fourier transformation of ^-functions in (|5.1| ) gives: 



t((Q)n,(P)n) = / n ^^ e^^^^PCCaik^aa)^), 
J k=i (2tt) 4 (2tt) 4 

where 

p((®l)N, {012) n) = «1,2---, ai,N\ «2,1, «2,2, Q!2,JV) 

has the form: 

p((ai)N, (« 2 )iv) = / fi{n ^%^e--^ x 

J ^=1 fc= i ra„! 

x 5 ^^e-^-}|a (n , m) | 2 . (5.2) 
k=i rnj 



Inserting ( |3."1| ) into (\b.'4 ) we find: 



p((a-)jy, (a + )jv) = exp{i E / dxda/[0 + (a;)Z} + _(:£ — 2/; a2,v)(f)-(x') — 

-0_(x)L>_ + (x - x'\ ai, i/ )0 + (x / )]}^(0+)^*(0-),(5.3) 

where 0_ is defined on the complex conjugate contour C_ : t —> t — is 
and D+( ), D -+(x — x'\ a) are the positive and negative frequency 

correlation functions correspondingly. 

We must integrate over sets (Q)n and (P)n if the distribution of fluxes 
momenta over cells is not fixed. In result, 

r(P) = J D A a 1 (P)d 4 a 2 (P)p((a 1 ) N , (a 2 ) N ), (5.4) 
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where the differential measure 

D 4 a(P) = ft P^K(P t (a) N ) 
takes into account the energy-momentum conservation laws: 

N N N 

K(P, (a) N ) = / J] d 4 Q„e l ^ a » Q »5^(P - £ Q 
The explicit integration gives that 

N 

K(P,(a) N )~ n<5 (3) («-«,), 



where a is the center of mass (CM) 3-vector. 

To simplify the consideration let us choose the CM frame and put a = 
(—i(3, 0). In result, 



poo N „ _ N 



K(E, (P) N ) = I™ n dE v e^^5{E - £ E, 

JU v=l v=l 

Correspondingly, in the CM frame, 

r(E) = f D/3 1 (E)Dp 2 (E)p((j3 1 )N i (J32)N) i 



where 

N dd 

Dm = n it:K(e, (j3) N ) 

and p((P)n) was defined in with = (— i/3fc i2 ,, 0), RePk,v > 0, k = 
1,2. 

We will calculate integrals over (3% using the stationary phase method. 
The equations for mostly probable values of (3 k : 

1 ° -K(E, (J3 k ) N ) = „ * xn f R((P)n), k = l,2, 



K(E,((3 k ) N )d(3 k> , v ' v ^ /iV/ R{(Pi)n) dfik,i 

(5.5) 

always have the unique positive solutions j3k )V (E). We propose that the 
fluctuations of j3 k near f3 k are small, i.e. are Gaussian. This is the basis of 
the local-equilibrium hypothesis [)32(| . In this case 1/ (3\ v is the temperature 
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in the initial state in the measurement cell v and l//?2,z/ is the temperature 
of the final state in the z/-th measurement cell. 

The last formulation QETgj) imply that the 4-momenta (Q)n and (-P)jv 
can not be measured. It is possible to consider another formulation also. 
For instance, we can suppose that the initial set (Q)n is fixed (measured) 
but (P) n is not. In this case we will have mixed experiment: (3\ iV is defined 
by the equation: 

E " = 'Jim:* 

and $2,v is defined by second equation in (|575Q . 

Considering limit N — > oo the dimension of cells tends to zero. In 
this case we are forced by quantum uncertainty principle to propose that 
the 4-momenta sets (Q) and (P) are not fixed. This formulation becomes 
pure thermodynamical: we must assume that {(3\) and (fy) are measurable 
quantities. For instance, we can fix {(3\) and try to find (/%) as the function 
of total energy E and the functional of (fli). In this case eqs.( |5T5|) become 
the functional equations. 

In the considered microcanonical description the finiteness of tempera- 
ture does not touch the quantization mechanism. Really, one can see from 
(|5.ij| ) that all thermodynamical information is confined in the operator ex- 
ponent 

e N(<pt<l>j) _ TT TT 

the expansion of which describes the environment, and the 'mechanical' 
perturbations are described by the amplitude Z(<p). This factorization was 
achieved by introduction of auxiliary field (j) and is independent from the 
choice of boundary conditions, i.e. from the choice of considered systems 
environment. 

5.2. Wigner functions formalism 

We will use the Wigner functions formalism in the Carrusers-Zachariasen 
formulation |25| . For sake of generality the m into n particles transition 
will be considered. This will allow to include into consideration the heavy 
ion-ion collisions. 

In the previous section the generating functional R((f3) n) was calculated 
by means of dividing the 'measuring device' (calorimeter) on the TV cells. 
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It was assumed that the dimension of device cells tends to zero (N — ► oc). 
Now we will specify the cells coordinates using the Wigner's description. 

Let us introduce the distribution function F n which defines the proba- 
bility to find n particles with definite momentum and with arbitrary coor- 
dinates. This probabilities (cross sections) are usually measured in particle 
physics. The corresponding Fourier-transformed generating functional can 
be deduced from (|5~3| ): 

F{z t (p + ) Ni {p-) N )= n ne /dj(9) ^ (9)e " %/(5) ^ (9) ^ (9) x 

xZ(0 + )Z*(0_). (5.6) 

The variation of F over z^(q) generates corresponding distribution func- 
tions. One can interpret z^(q) as the local activity: the logarithm of z^(q) 
is conjugate to the particles number in the cell v with momentum q for the 
initial (ij = 21) or final (ij = 12) states. Note that z"Aq)4>*(q)(j)j(q) can 
be considered as the operator of activity. 

The Boltzman factor e~ l3i ^ e ^ can be interpreted as the probability to 
find a particle with the energy e(q) in the final state (i = 2) and in the 
initial state (i = 1). The total probability, i.e. the process of creation 
and further absorption of n particles, is defined by multiplication of this 
factors. 

The generating functional ( |5.tjj ) is normalized as follows: 



F(z = l,((3)) = R(((3)), (5.7) 
F(z = 1 (f3)) = \Z(O)\ 2 = po(0 ± )U ±= o 

Where 

po((f>±) = Z(<j> + )Z*(<J>-) 

is the 'probability' of the vacuum into vacuum transition in presence of 
auxiliary fields <f>±. The one-particle distribution function 

Fi((MN,(P2) N ;q) = jj^F\ z=0 = 

= {0*(g)e-^^/ 2 }{0 J (g)e-^^/ 2 } /?o (0 ± ) (5.8) 
describes the probability to find one particle in the vacuum. 
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So, 

= /dr^e^e^^ (5-9) 



We introduce using this definition the one-particle Wigner function W\ 



^i((A)iV, (Mn; r,q) = J dye^e-^Ur + y/2)^(r - y/2)p (<f> ± ). 



This distribution function describes the probability to find in the vacuum 
particle with momentum q at the point r in the cell v 

Since the choice of the device coordinates is in our hands it is natural 
to adjust the cell coordinate to the coordinate of measurement r: 



W 1 ((p 1 ) Nl (p 2 ) N ; r,q) = J dyeWe-M'MyUr + y/2)^(r - y/2)p (cf> ± ). 



This choice of the device coordinates lead to the following generating func- 
tional: 



F(z, 0) = exp{i J dydr[4>+(r + y/2)D+_(y; (3 2 (r), z)$-(r - y/2) - 

-0_(r + y/2)D- + (y; A(r), z)4>+(r - y/2)] }p (^±), (5.10) 



are the modified positive and negative correlation functions. 

The inclusive, partial, distribution functions are familiar in the particle 
physics. This functions describe the distributions in presence of arbitrary 
number of other particles. For instance, one-particle partial distribution 
function 



Fi((Pi)n, (fo) N \ q) == J drW^N, ((3 2 ) N ; r, q). 



So 



where 




Pij(r,q; (/?)) 



S 



5zij(r, q) 



z=l — 



e -0i(r)e(q) 



I 



dye^Ur + y/2)^(r - y/2)p(<f> ± , (/?)) 



(5.11) 



(27T)3 e (g) 
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where eq.(|5T7|) was used. 

The mean multiplicity n^-(r, q) of particles in the infinitesimal cell Y 
with momentum q is 

n ij (r,q) = J dq ^ ^ In F(z, (f3))\ z=i . 

If the interactions among fields are switched out we can find that (omitting 
indexes): 

n(Y,q ) = e/3(r)qo _ r q = e(q)>0. 

This is the mean multiplicity of black-body radiation. 
5.3. Closed path boundary conditions. 

The developed formalism allows to introduce more general 'closed-path' 
boundary conditions. Presence of external black-body radiation flow will 
reorganize the differential operator exp{iV(0*0j)} only and new generating 
functional p cp has the form: 

Pc P (a 1 ,a 2 )=e*WMpo(<l>±). 

The calculation of operator N((p*(f)j) is strictly the same as in Sec. 3. Intro- 
ducing the cells we will find that 

Nfflfa) = J drdyUr + y /2)h l] {Y,y)^{r - y/2), 

where the occupation number hij carries the cells index r: 

hij(r, y) = J duj(q)e iqy n tj (r, q) 

and (q = e(q)) 

n ++ (r,q ) = n__(r,g ) = n(r, (A + f3 2 )\q \/2) = e[ ^ m \ qa \, 2 _ r 

n+-(r,q ) = Q(q )(l + n(r,(3 2 q )) + Q(-qo)n(r, -faqo), 
n-+(r, q ) = n+_(r, -q ). 
For simplicity the CM system was used. 
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Calculating po perturbatively we will find that 



Pcp(P) = exp{-iV(-ij+) + iV(-ij-)} x 



exp{i J drdy[ji(r + y/2)G t] (y, (P(r))jj(r - y/2)} 



(5.12) 



where, using the matrix notations, 




iG(q, (J3(r))) = « 2 - 




+ 



) 



(5.13) 



and 



= _ e /3(ko|±go)/2_ 



(5.14) 



Formally this Green functions obey the standard equations in the y space: 



since &((7oo) ^ reflects the mass-shell particles. But the boundary con- 
ditions for this equations are not evident. 

It should be underlined that in our consideration r is the coordinate of 
measurement, i.e. r is as the calorimeter cells coordinate and there is not 
necessity to divide the interaction region of QGP on domains (cells). This 
means that L must be smaller then the typical range of fluctuations of 
QGP. But, on other hand, L can not be arbitrary small since this will lead 
to assumption of local factorization property of correlators, i.e. to absence 
of interactions. 

So, changing (3 — > (3{r) we should assume that A(2)(0 an d (r, k) 

are constants on interval L. This prescription adopts Wigner functions 
formalism for the case of high multiplicities. It describes the temperature 
fluctuations larger then L and averages the fluctuations smaller then L 
leading to absence, in average, of 'non- Gaussian' fluctuations. 

It is the typical 'calorimetric' measurement since in a dominant number 
of calorimeter cells the measured mean values of energy, with exponential 
accuracy, are the 'good' parameters ~ 1/^2 (r, E). We will assume that 



{d 2 - m 2 ) y G lt = %), 
(d 2 - m^yGij = 0, i^j 
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the dimension of calorimeter cells L « L cr , where L cr is the dimension 
of characteristic fluctuations at given n. In deep asymptotic over n we 
must have L cr — ► oo. This consideration shows that the offered experiment 
with calorimeter as the measuring device of particles energies is sufficiently 
informative in the high multiplicities domain. 



6 Nonstationary statistical operator 

One can not expect the evident connection between the above considered 
S'-matrix (microcanonical) and Zubarev's [j32j approaches. The reason is 
introduction into Zubarev's formalism interaction with a heat bath, exter- 
nal to system under consideration. This interaction is crucial for definition 
of NSL for explanation of the trend to maximal-entropy state, starting 
evolution from local-equilibrium stateQ. 

Therefore, in Zubarev's theory the local-equilibrium state was chosen as 
the boundary condition. It is assumed that in the suitably defined cells of 
the system at a given temperature distribution T(x,t) = l/(3(x,t), where 
(x, t) is the index of the cell, the entropy is maximal. The corresponding 
nonequilibrium statistical operator 

describes evolution of a system in the time scale t. Here is the energy- 
momentum tensor. It is assumed that the system 'follows' to (3(x, t) evolu- 
tion and the local temperature T(x, t) is defined as the external parameter 
which is the regulator of systems dynamics. For this purpose the special 
ze-prescription was introduced (it was not shown in (|6.1|) [j52fl . It brings 
the interaction with heat bath. 

The KMS periodic boundary condition can not be applied for nonsta- 
tionary temperature distribution and by this reason the decomposition: 

P(x,t) =p + (3 1 (x,t) (6.2) 

was offered in the paper [|3D|. Here /3q is the constant and the inequality 

A) » \Pi(x,t)\ 

7 This condition is not necessary in the ^-matrix formalism since it is 'dynamical' by its nature, i.e. 
includes the notion of initial- and final-states as the boundary conditions. 
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is assumed. Then, 

pz „ e -MH 0+ V + B) 3) 

where Hq is the free part of the Hamiltonian, V describes the interactions 
and the linear over (3i/(3q term B is connected with the deviation of tem- 
perature from the 'equilibrium' value l//?o- Presence of 5-perturbations 
creates the 'thermal' flows in the system to explain increasing entropy. 
Considering V and B as the perturbations one can calculate the observ- 
ables averaging over equilibrium states, i.e. adopting the KMS boundary 
condition. Using standard terminology one can consider V as the 'mechan- 
ical' and B as the 'thermal' perturbations. 

The quantization problem of operator ( |6.3|) is connected with definition 
of the space-time sequence of mechanical (V) and thermal (B) excitations. 
It is necessary since the mechanical excitations give the influence on the 
thermal ones and vice versa. It was assumed in that V and B are 



commuting operators, i.e. the sequence of V- and B-perturbations is not 
sufficient. The corresponding generating functional has the form 

Z(j) = exp{-i f d A x{V{-ij{x)) + ^^-T 00 [-ij(x 

[Dp 



oo In 



A) 

where the time contour Cp was described in Sec. 4.1 and r is the measure- 
ment time. 

It is evident that this solution leads to the renormalization by the in- 
teractions with the external field j3(x, t) even without interactions among 
fundamental fields The source of this renormalizations is the kinetic 
term in the energy-momentum tensor Too, i- e - follows from 'thermal' in- 
teractions with external heat bath. Note absence of this renormalizations 
in the ^-matrix formalism, see, for instance (|5.2(j|) , where the interactions 
generated by V-perturbations only. 



In |53j] the operators V and B are noncommuting ones and 5-perturbations 



were switched on after ^-perturbations. In this formulation the nondy- 
namical renormalization are also present but it is not unlikely that they 
are canceled at the very end of calculations [|53 . 
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This formulation with (3(x, t) as the external field reminds the old, firstly 
quantized, field theory in which matter is quantized but fields are not. It 
is well known that consistent quantum field theory requires the second 
quantization. Following to this analogy, if we want to take into account 
consistently the reciprocal influence of V- and 5-perturbations the field 
(3(x,t) must be fundamental, i.e. must be quantized (and the assumption 
of paper []3D[ becomes true). But it is evidently the wrong idea in the 



canonical Gibbs formalism. So, as in the firstly quantized theory, the 
theory with operator must have the restricted range of validity 



7 Conclusion 

In our interpretation of the real-time finite-temperature field theory the 
statistics and the fields quantum dynamics were factorized: statistics is 
fixed by the operator exp{iV(0*0j} and a pure field-theoretical dynamics 
is described by po((f>±) = Z(0 + )Z*(0_), where Z((f)±) is the vacuum into 
vacuum transition amplitude in presence of the external (auxiliary) fields 
< vac\vac ></,. We can say that the operator exp{7V(0*0j} maps the system 
of interacting fields on the state with definite thermodynamical parameters. 
We had concentrated our attention in this paper on the structure and 
origin of operator exp{N((j)*(j)j} only and do not discuss po(4>±). But the 
developed formalism allows to use following '^-matrix' properties which 
are new for thermodynamics to define Pq. 

First of them is absence of Matsubara imaginary parts of time contour in 
Pq by definition: the approach is pure 'real-time'. This allows to construct 
the formalism without referring to time asymptotic properties of correla- 
tion (Green) functions, and introduce the temperature description without 
using a notion of grand canonical ensemble constructing the environment 
of the system, i.e. the measuring device, 'by hand'. 

Moreover, discussed factorization property have important consequence 
which would allow to calculate expectation values with high accuracy. Let 
us consider the theoretical problem of the po(4>±) calculation. To define the 
functional measure the ortho-normalizability (i.e. the unitarity) condition 
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may be used. It leads to following representation 

p o (0) = e~ ik ^ J DM(<$>)e~ u ^ e) ef dx{v '^ +j) t (7.1] 
where the expansion over operator 

5 5 



5j(x) 5e(x) 



K(j, e) = 2Re J dx 

generates perturbation series and 

e ) = V($ + e) - V($ - e) - 2Re J dxev'(§) 

weights quantum fluctuations. The most important term in (|Y.1|) is the 
measure 

M($) = n d$(x)5(d 2 J> + m 2 $ + v'{<!>) - j) 

X 

where v'(§) = 6V(Q>)/5$>(x). So, solving the equation 

#2$ + m 2 $ + ?/($) = j (7.2) 

we will find all contributions^. 

At the very end of calculations one must put e = j = 0. Therefore, 
eq . ( [7T2D can be solved expending it over j. This shows that ( |7TT|) restores 
at j = the usual stationary phase method. Indeed, it can be verified that 
Qy.-U ) gives usual perturbation theory . 

But the eq.( frr^) gives much more possibilities. Note that l.h.s. of this 
equation is sum of known classical forces and the r.h.s. is the quantum 
force j. Eq. (\(.'4) establish the local equilibrium between this forces. This 
solves the old standing problem of quantization with constraints: it can be 
done by field transformations in path integrals since the eq. (|772]) shows the 



way as j must be transformed when the l.h.s. is transformed. Presence of 
derivatives in ( \(.'4) shows that the quantum force must be transformed in 



the tangent space of fieldsQ. 

The r.h.s. of eq. (\(.'4 ) may contain also an additional force to describe 
the external influence on the system of interacting fields. This force was 

8 This means that the unitarity condition is necessary and sufficient for definition of path integral 
measure for Ro((p±) |56 

9 This explains why the ordinary transformation of path integral is impossible, gives wrong result [b7j. 
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omitted in eq. flTTg ) assuming that a process of particles creation (and ab- 
sorption) switched on adiabatically. 

As was mentioned above the action of operator e~ N ^ ,z ' ,c ^> on Ro(4>) maps 
interacting fields system on measurable states. Let us consider what this 
gives. Result of action has the form: 

p(J3,z) = e ~ ik ^ J DM(<S>)e- u ^e- N ^ z -^\ 

where N = N\ + N2 and 

N m {J3,z;$) = / «a;(A:)e^)W^2 + _ ( _ +) (A:,r)|r(A:,$)| 2 . (7.3) 

Here r is considered as the index of calorimeter cell. This formulae needs 
more careful explanation. Instead of (|Y.3| ) we must consider 

J dq5 L (q)T(k + q, <S>)T*(k - q, $). (7.4) 

where L is the scale where (3i(o\(r) and z + _t_ + \(k,r) can be considered as 
the constants (L is the dimension of calorimeter cell). If L — > 00 then 5i{q) 
can be changed on usual ^-function 5(q) and, therefore, in this limit we will 
have Q/.4 ). We had considered this limit supposing that the measurement 



is not in contradiction with quantum uncertainty principle. 

So, deriving Nu2)(P, z; $) there was used the condition that r is the 
coordinate of size L cell. With this condition 

T(k, $) = J dxe lkx {dl + m 2 )$ (7.5) 

can be considered as the order parameter. Indeed, V(k, $) is the element 
of actions symmetry group since it is linear over field $ and the generating 
functional R((3, z) is trivial if < \T(k, $)| 2 >= 0. In this case there is not 
creation of particles, i.e. there is not measurable asymptotic states (fields). 
Indeed, it can bee shown [|58j] that all quantum corrections to solitons 



contribution in the (l+l)-dimensional sin-Gordon model equal to zero. 
This is in accordance with result of |59j and with factorizability of solitons 



(S-matrix. Then it is easily seen computing integral in ( |7.5|) by parts that 
T(k, $ s ) = 0, where $ s is the soliton solution. This result shows that hid- 
den symmetry of sin-Gordon model can not be broken and corresponding 
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(polynomial) integrals of motion are conserved. The application of this 
idea for non-Abelian field theory should be fruitful. 
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